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Summary. — 

We show analytically that GRB jets with an emissivity profile given by the Fisher 
distribution, e(6) = A ■ exp(B ■ cos 8), have the unique property of producing equal 
numbers of bursts per logarithmic interval in Ei so , and therefore in most burst 
properties. Since this broad distribution of burst properties is a key feature found 
by HETE-2, a Fisher-shaped universal jet model can explain many of the observed 
properties of XRFs, X-ray-rich GRBs, and GRBs reasonably well, in contrast to 
a power-law universal model. For small viewing angles, the Fisher distribution 
can be approximated by a Gaussian, whose properties have been explored by |12) . 
We also show that the Fisher universal jet model produces a broad distribution in 
the inferred radiated energy -E^ nf , in contrast to the narrow distribution predicted 
by the uniform variable opening-angle jet model UJ. Here we present Monte Carlo 
simulations of both a Fisher-shaped universal jet model and a Fisher-shaped variable 
opening-angle jet model. 

PACS 98.70.Rz - gamma-ray sources; gamma-ray bursts. 



1. Introduction 



The HETE-2 results show that the properties of XRFs [Sj, X-ray-rich GRBs, and 
GRBs form a continuum in the [S E (2 - 400 kev), E°^l k ] plane [IB]. They also show that 
the relation between the isotropic-equivalent burst energy E- lso and the peak energy -E pea k 
of the burst spectrum in vF v in the rest frame of the burst found by extends to XRFs 
and X-ray-rich GRBs [5] . A key feature of the distribution of bursts in these two planes 
is that the density of bursts is roughly constant along these relations, implying equal 
numbers of bursts per logarithmic interval in Se, -Speak) -^iso and -E pea k- These results, 
when combined with earlier results (SHE], strongly suggest that all three kinds of bursts 
are the same phenomenon. It is this possibility that motivates us to seek a unified jet 
model of XRFs, X-ray-rich GRBs, and GRBs. 

In our previous paper [5], we explored two different phenomenological jet models: a 
variable jet opening-angle model in which the emissivity is uniform across the surface 
of the jet and a universal jet model in which the emissivity is a power-law function of 
the angle relative to the jet axis. We showed that while the variable jet opening-angle 
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model can account for the observed properties of all three kinds of bursts, the power-law 
universal jet model cannot easily be extended to account for the observed properties 
of both XRFs and GRBs. In response to that conclusion, El considered a quasi- 
universal Gaussian jet model They showed that such a model can explain many of 
the observed properties of XRFs, X-ray-rich GRBs, and GRBs reasonably well. 

Here we consider a universal jet model in which the emissivity of the jet as a function 
of viewing angle is a Fisher distribution (such a distribution is the natural extension 
of the Gaussian distribution to a sphere). We show that the Fisher distribution has 
the unique property that it produces equal numbers of bursts per logarithmic interval in 
.Eiso, and therefore in most burst properties, consistent with the HETE-2 results. We also 
show that the Fisher universal jet model produces a broad distribution in the inferred 
radiated energy E^. This is not the case for variable opening-angle jet models because 
of the extra degree of freedom provided by the distribution of opening angles. Thus we 
find that the Fisher universal jet model considered here and the variable opening-angle 
jet model discussed in [§] make different predictions for the distribution in Eif. Further 
observations of XRFs can determine this distribution and therefore distinguish between 
these two models of jet structure. For completeness, we also simulate a variable opening- 
angle jet model whose emissivity profile is also a Fisher distribution, and we find similar 
results to the universal Fisher jet model. 

2. Simulations 

The Fisher distribution is the only universal jet profile that satisfies the following 
two constraints. Let X- lso = ln_E; so , and fi v = cos# v , where 9 V is the angle between the 
line-of-sight and the center of the jet. Observations tell us that roughly dN/dX iso = oi, 
where a\ is some constant. By the definition of 8 V we know that dN/dfi v = a,2, where 
ai is another constant. We can describe a universal jet as an arbitrary function of /z v : 
^iso = /(a*v)- We wish to choose f(fi v ) so as to satisfy these constraints. 

Note that 

dN dN dX iso 
1 ° 2 = T~ = IV 1 ' 

and therefore df/dfj, v = a 2 /&i, integrating this expression gives the Fisher distribution 

(2) E iso = 4ttA ■ e B - cos0 v = 4ttA ■ e (^^-i)/d 2 ^ 

where B = 6$ 2 . In the small 9 V limit, this reduces to a Gaussian jet, Ei so = Air A -e^/ 2e o. 
Integrating the emissivity over the jet gives E^ rne , 

f 71 / 2 4-wA 

(3) £* ruc = 2 • 2nA / e B - cose sinOdO = (e B - 1). 

In this work we consider two oppositely directed jets, hence the leading factor of 2. 
We note that for any non-uniform jet this quantity is not the same as the E 1 inferred 
using the method outlined by That quantity we term, E" 11 = Ei so ■ (1 — cos^j), 
where } = max(6*o,0 v ) j210- For #o = 0.1 rad, the quantities E- mo and E™ 1 vary over 
domains of ~ 43 and ~ 40 decades, respectively, although observational selection effects 
will truncate both of these distributions. We perform Monte Carlo simulations using the 
method presented in using the detector thresholds from the WXM on HETE-2. 
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3. Results 

We consider 3 models. (1) A Universal Fisher Jet model with log.E* ruc = 51.1 and 8q 
values drawn from a log-normal distribution with width 0.2 and log#g — —1-0 (UFJ1), 
following the parameters of 0^], (2) a Universal Fisher Jet model with log#Q = —1.3 
and log£* rue = 51.8 (UFJ2), and (3) a Variable Opening- Angle Fisher (VOAFJ) jet 
model with log-E^ ruc = 51.5 and 9q values drawn from a power-law with index apL = 
—3.3 and extending for two decades from a maximum of tt/2. 

Figures 1 and 2 show the distribution of detected and non-detected bursts in various 
planes for our 3 models. All 3 models exhibit roughly equal numbers per logarithmic 
decade in E lso . For the universal models, this is a natural consequence of the Fisher profile 
of the jet. For the VOAFJ model this is a consequence of our choice of apL = —3.3. We 
note that the UFJ1 model is unable to accomodate the highest observed values of Ei so and 
Epcak- The maximum Ei so generated by a Fisher jet is approximately E™ x = E^ uc /0q, 
which gives - 1.2 x 10 53 erg for UFJ1 and - 2.5 x 10 54 erg for UFJ2. 




Fig. 1. - Scatter plots of detected (black) and undetected bursts (gray) in the [9o,9 v ] plane (top 
row) and [E- lso ,9 v ] plane (bottom row) for our 3 models: UFJ1 (left), UFJ2 (middle) and VOAFJ 
(right). The triangular region in the upper-left corner represents bursts that we do not simulate 
to increase the percentage of detected bursts in a sample of 50,000. 

The bottom row of Figure 2 shows the histogram of E™ 1 values for the detected 
bursts. It is clear that the E™ distribution does not agree with the inputted J5'^ rue 
distribution. For example, in the UFJ1 model, the peak of the E^ ue distribution was 
chosen to correspond to the "standard energy" found by |2] , however the model is 
unable to recover that value using their method. The E" 11 distribution typically peaks 
at a lower energy and has a tail extending to even lower energies. It may be the case 
that the observed distribution of E™ { values does extend down to lower energies. 
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Fig. 2. - Top row shows scatter plots of detected (black) and undetected bursts (gray) in the 
[£iso,-E pcak ] plane for our 3 models: UFJ1 (left), UFJ2 (middle) and VOAFJ (right). Bottom 
row shows histogram of E 1 ™ 1 values for detected bursts, as compared to the input distribution 
of E\ ruc (dashed curve) . 



4. Conclusions 

Both universal and variable Fisher jet models can be found that reproduce most of 
the observed properties of XRFs and GRBs. To accomodate the highest -E^o bursts, very 
small jet opening-angles (~ 2° — 3°) may be required in both the variable opening-angle 
uniform jet models and in the Fisher models. Finally, may be a powerful probe 
of jet structure, as various models give different predictions for its distribution. More 
observations of XRFs with redshifts and jet-break times are crucial to answering this 
question, highlighting the importance of continuing HETE-2 during the Swift mission. 



REFERENCES 

[I] Amati, L., et al., A&A, 390 (2002) 81. 

[2] Bloom, J., Frail, D. A., & Kulkarni, S. R., ApJ, 588 (2003) 945. 
[3] Dai, X. & Zhang, B., ApJ, submitted, 2004 (astro-ph/0407272 ). 
[4] Frail, D., et al., ApJ, 562 (2001) L55. 

[5] Heise, J., et AL., 2000, in Proc. 2nd Rome Workshop: Gamma-Ray Bursts in the Afterglow 
Era, eds. E. COSTA, F. Frontera, J. Hjorth (Berlin: Springer- Verlag) , 16. 

[6] Kippen, R. M., et AL., 2002, in Proc. Gamma-Ray Burst and Afterglow Astronomy, AIP 
Conf. Proceedings 662, eds. G. R. Ricker & R. K. Vanderspek (New York: AIP), 244. 

[7] Kumar, P. & Granot, J., ApJ, 591 (2003) 1075. 

[8] Lamb, D. Q., et al., ApJ, submitted, 2005. 

[9] Lamb, D. Q., Donaghy, T. Q., and Graziani, C, ApJ, 620 (2005) 355. 
[10] Sakamoto, T., et al., ApJ, submitted, 2004b (astro-ph/0409128 ). 

[II] Zhang, B., & Meszaros, P., ApJ, 571 (2002) 876. 

[12] Zhang, B., Dai, X., Lloyd-Ronning, N. M., & Meszaros, P., ApJ, 601 (2004) L119. 



